Full-Coupled Channel Approach to Doubly Strange s-Shell Hypernuclei 
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We describe ab initio calculations of doubly strange, S = —2, s-shell hypernuclei (aaH, aaH, Aj\H e 
and AAHe) as a first attempt to explore the few-body problem of the /wZZ-coupled channel scheme 
for these systems. The wave function includes AA, AE, NS and EE channels. Minnesota AW, D2' 
YN, and simulated Y Y potentials based on the Nijmegen hard-core model, are used. Bound-state 
solutions of these systems are obtained. We find that a set of phenomenological B$Bg interactions 
among the octet baryons in S = 0, — 1 and —2 sectors, which is consistent with all of the available 
experimental binding energies of S = 0, — 1 and —2 s-shell (hyper-)nuclei, can predict a particle 
stable bound-state of aaH. For aaH and AAHe> AA — EAT and SN — AE potentials enhance the net 
AA — ATE coupling, and a large 3 probability is obtained even for a weaker AA — AE potential. 

PACS numbers: 21.80.+a, 21.45,+v, 21.10.Dr, 13.75.Ev 



Both recent experimental and theoretical studies of 
doubly strange (S = —2) s-shell hypernuclei ( AA H, AA H, 
AA He, and AA He) are the utmost exciting topics in the 
field of hypernuclei HSHSHSli[lfl[l3 An 
experimental report lj on a new observation of AA He has 
had a significant impact on strangeness nuclear physics. 
The Nagara event provides unambiguous identification 
of AA He production, and suggests that the AA interac- 
tion strength is rather weaker than that expected from 
an older experiment |l2j. 

The BNL-AGS E906 experiment has conjectured 
a formation of » f H, in accordance with our earlier 
predictions [l3l|l4j that AA H would exist as a particle sta- 
ble bound-state against strong decay. If this is the case, 
the AA H would be the lightest bound state among doubly 
strange hypernuclei. However, a theoretical study|3j of 
the weak-decay modes from AA H does not support this 
conjecture, and our earlier studies should be reanalyzed 
by taking account of the new datum, Nagara event. 

A recent Faddeev-Yakubovsky search for AA H0 found 
no bound-state solution over a wide range of AA inter- 
action strengths, although this conclusion has been in 
conflict with the result calculated by authors using a vari- 
ational method^. The total binding energy is more sen- 
sitive to the 3 Si channel of the AN interaction than to 
the 1 So AA interaction, because the number of the 3 Si 
AN pairs is three times larger than the number of the 
1 Sq AA pair, as was discussed in Ref. Q. Therefore, the 
spin-dependent part of the AN interaction has to be de- 
termined very carefully. The algebraic structure of the 
(<r A -(Tn) interaction for the S = —2 systems is similar to 
the structure for the A He. Namely, the AA" interaction, 
which is utilized in the theoretical search for AA H, has to 
reproduce the experimental B A ( A He) as well as the B\s 
of A = 3, 4 S = — 1 hypernuclei. However, there is a long- 
standing problem known as the A He anomaly 15J, since 



the publication by Dalitz et al. in 19720. Recently, 
Akaishi et «L|17l| successfully resolved the anomaly by 
explicitly taking account of AA^ — EA^ coupling. 

Considering the fact that the AA system couples to 
the NE and EE states, and also the AN system couples 
to the EA^ states, a theoretical search for AA H should 
be made in a fully coupled channel formulation with a 
set of interactions among the octet baryons. The AA H- 
IH (or AA He-|He) mixing due to AA— NE coupling is 
also interesting topic, since the a-formation effect could 
be significant pi li| Thus, the purpose of this study is 
threefold: First is to describe a systematic study for the 
complete set of s-shell hypernuclei with S=—2 in a frame- 
work of a full-coupled channel formulation. Second is to 
make a conclusion if a set of baryon-baryon interactions, 
which is consistent with the experimental data, predicts 
a particle stable bound state of AA H. The third is to ex- 
plore the fully hyperonic mixing of AA H, including the 
AN— UN transition potential in addition to AA— NS— EE. 

The wave function of a system with S — —2, compris- 
ing A(— N + Y) octet baryons, has four isospin-basis 
components. For example, AA He has four components as 
ppnnAA, NNNNNE, NNNNAT, and AW AWES. We 
abbreviate these components as AA, A^S, AE and EE, 
referring to the last two baryons. The Hamiltonian of 
the system is hence given by 4 x 4 components as 
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where Hb 1 b 2 operates on the -B1-B2 component, and 
Vb 1 b->-b[ b' 2 is the sum of all possible two-body transition 
potential connecting the B\Bi and B' 1 B' 2 components: 
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Va.A-AS = y^^ATjA-AriS; Vjvh-AS = Ujvh-AE> (3) 
i=l 

VAA-ES = WAA-SS, V/VE-EE = Wath-SE, (4) 

VaE-EE = ^ VNjA-NiX + ^AE-SS- (5) 
»=1 

Note that we take account of /^/-coupled channel po- 
tentials including the AN— UN and ATH— AE transitions 
(Eq.||3J|) in the 3 Si channel , wh ile other full-coupled chan- 
nel approaches (e.g., Refs. [lollll| ) only take AA-VS-EE 
in the ^So channel (Eqs.(J2J and J2J) into account. 

In the present calculations, we use the Minnesota 
potential|l£| for the NN interaction and D2' for the 
YN interaction. The Minnesota potential reproduces 
reasonably well both the binding energies and sizes of 
few-nucleon systems, such as 2 H, 3 H, 3 He and 4 He[l9j. 
The D2' potential is a modified potential from the orig- 
inal D2 potential^?!- The strength of the long-range 
part (V b in Table I of Ref. of the D2' potential in 
the AN- AN 3 Si channel is reduced by multiplying by 
a factor (0.954) in order to reproduce the experimental 
BA{%He) value. The calculated Ba values for the A hy- 
pernuclei ( A H, a H, A He, A H*, A He*, and A He) are 0.056, 
2.23, 2.17, 0.91, 0.89, and 3.18 MeV, respectively. For the 
YY interaction, we use a full-coupled channel potential 
among the octet baryons in both the spin triplet and the 
spin singlet channels. We assume that the YY potential 
consists of only the central component, and the effect 
due to the non-central force (e.g., tensor force) should 
be included into the central part effectively. The YY 
potential has Gaussian form factors, whose parameters 
are set to reproduce the low-energy S matrix of the Ni- 
jmegen hard-core model D (ND) or F (NF)|20T|. We take 
the hard-core radius to be r c = 0.56271(0.44915) fm in 
the spin singlet (triplet) channel for the ND, whereas 
r c = 0.52972(0.52433) fm is used in the singlet (triplet) 
channel for the NF. Each number is the same as the hard- 
core radius of the YN sector in each channel for each 
model. The strength parameters are firstly determined 
on a charge basis, and then the strength parameters on 
an isospin basis are constructed from the charge-basis pa- 
rameters. We denote NDs (NFs) for the simulating ND 
(NF) potential. 

The calculations are made by using the stochastic vari- 
ational method|2lll22j|. This is essentially along the lines 
of Ref. [23|], except for the isospin function. The isospin 
function consists of four components, in accordance with 
Eq. (JTJ. The reader is referred to Refs. jU for the 
details of the method. 

Table 0] lists the Baa values for S = —2 hypernuclei. 
Using NDs or NFs YY potential, we have obtained the 
bound-state solutions of AA H, AA H, AA He and AA He. In 
the case of ND5 YY potential, we have 

AB A C A alc) UHe) = i?i C A alc )( A lHe)-2i?i Calc )(iHe) 



TABLE I: AA separation energies, given in units of MeV, of 
A = 4 — 6, S — —2 s-shell hypernuclei. 
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Baa(aIH) Baa(aIH) 


BAA( A lHe) 


BAA( A |He) 


ND S 


0.107 4.04 


3.96 


7.93 


niND s 


0.058 3.74 


3.66 


7.53 


NF S 


0.127 3.84 


3.77 


7.52 


Exp 






7.25 ± 0.191^? 



= 1.58MeV, (6) 

which is slightly larger than the experimental value 

(A£ AA xp )( A ^He) = 1.01 ± 0.20±°;^ MeV lj), while the 
NFg reproduces the experimental value fairly well. We 
have also calculated these systems using a modified NDs 
(mND s ) YY potential; The strength of the A A diago- 
nal part of the mNDs potential is reduced by multiply- 
ing by a factor of 0.8 from the original NDs in order to 
reproduce the experimental AB\\( A ^He). The scatter- 
ing length and effective range parameters for the NDs, 
mNDs, and NFs are (given in units of fm): (a s ,r s ) = 
(-1.37,4.98), (-0.91,6.25), and (-0.40,12.13), respec- 
tively. The scattering length for the mNDs or NFs is 
consistent with the other analyses 0, ^} concerning the 
Nagara event. We should note that the mNDs potential 
predicts the particle stable bound state of AA H; The ob- 
tained energy is very close to, but still (0.02 MeV) lower 
than, the A H + A threshold. Therefore, due to the result 
for mNDs or NFs > we should come to the following novel 
conclusion: A set of phenomenological baryon-baryon in- 
teractions among the octet baryons in S = 0,-1 and 
—2 sectors, which is consistent with the Nagara event as 
well as all the experimental binding energies of S — and 

— 1 s-shcll (hyper-) nuclei, can predict a particle stable 
bound state of AA H. 

Figure H schematically displays the present results of 
the full-coupled channel calculations oi A = 3 — 6,5 = 

— 1, —2 hypernuclei, using the mNDs- Since the present 
calculation has been made on the isospin basis, the re- 
sults for AA He are qualitatively similar to the results 
for AA H, so that we omit the explicit result for AA He. 
Fig. H also displays the probabilities of the NS, AE 
and EE components for the S = —2 hypernuclei. In 
the case of NFs, the probabilities are (given in per- 
centage): (Pjvh,Pae,Pse) = (0.58,0.37,0.03) for AA H, 
(3.10, 2.10, 0.10) for ^H, and (1.33, 1.13, 0.10) for AA He, 
respectively. In the present calculations, the AA compo- 
nent is the main part of the wave function. No unrealistic 
bound states were found for the YY subsystem, since the 
hard-core model hardly incorporates an unrealistic strong 
attractive force in the short-range region, in contrast to 
the soft core model, such as NSC97e or NSC97f[ll|. This 
is one of the reasons why we used the YY potential con- 
structed from the hard-core model, for the first attempt 
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FIG. 1: A and AA separation energies of A = 3 — 6, S = — 1 and -2 s-shell hypernuclei. The Minnesota NN, D2' YiV and 
mNDs YY potentials are used. The width of the line for the experimental Ba or Baa value indicates the experimental errorbar. 
The probabilities of the N3, AE and EE components are also shown for the AA hypernuclei. 



to the full-coupled channel calculation. Pss's are very 
small for both systems, due to a large mass difference be- 
tween the AA and EE channels (toes — jraa — 155 MeV). 

We should emphasize that the Pns(aa^) obtained 
by the mNDs has a surprisingly large value (4.55%), 
which is larger than the Pjvh(aaH) obtained by the NFs 
(3.10%), in spite of the fact that the strength of the 
AA — NE coupling potential of the ND is rather weaker 
than that of the NF. This does not imply that a stronger 
AA— NE coupling potential means a larger Pns probabil- 
ity. This is in remarkable contrast with other calculations 
based on (i+A+A) and (a+E~) two-channel model0,@- 

Although the present calculation assumes no simpli- 
fied structures, such as (t + A + A) and (a + H~), this 
kind of model is useful to make a clear explanation 
of the complicated full coupling dynamics of the A — 
5,3=— 2 hypernucleus. Let us consider a set of simple 
core nucleus +Y(+Y) model wave functions for the 
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| A ^H) =1p t X l/'AA X VAA-t, 
||H) = lp a X lp s - X ^ H -_ a , 

T 



(7) 
(8) 



Ia! H ).Sas = \/£ VPt X [V'ASoIsaJ X 



3 bPh X [^AS-]s AS ] X 1pAH--h 

(for 5 AS = or 1), (9) 



where ip c (c = t, h, a) is the wave function (WF) of the 
core nucleus, ifj Y Y (YY = AA, 2~,AE) is the WF of 
the hyperon(s), and ipYY-c is the WF that describes the 



relative motion between YY and c. We assume that all 
of the baryons occupy the same (0s) orbit. For the ae-H 
state, we have two independent states for the WF V'as, 
that the spin of two hyperons (Sad) is either a singlet or 
a triplet. Since the EE component plays a minor role, we 
omit the E |.H state. Using these WFs, we can obtain the 
algebraic factors for each averaged coupling potential of 
the allowed spin state, v 3 or v l : 



(Vaa-ate) = 
(Vaa-ae) = 

(Vns-at.) = • 
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(for 5 A2 = 1), 



!^a~-ae ( fOT Sax = 0), 
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The vaa-ne potential is suppressed by a factor of yl/2 
for the A — 5 hypernucleus. The vna-ni: and «ae-ae 
potentials, particularly in the spin triplet channel, play 
significant roles instead. Namely, these equations imply 
that the AE component strongly couples both to the AA 
and to the NE components, and the AE component plays 
a crucial role in the hypernucleus. 

The normalized energy expectation values of the 
Hamiltonian for AA \H are (given in units of MeV), 
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(for the mNDc 



(13) 



(for the NF 5 ). 



\ -14.92 -10.01 



Here, we display only the 3x3 components of the Hamil- 
tonian 0J, comprising AA, NE, and AS, since the con- 
tributions from the EE component are not large. If wc 
solve the eigenvalue problem, det(h — XI) = 0, we obtain 
the ground state energy, £ = -11.82 MeV (-11.82 MeV), 
and the probability, P NE = 3.99% (2.83%), for the mND 5 
(NFg). The first 2x2 components for the mNDj are quite 
different from those for the NFs, while the last row and 
the last column are qualitatively similar to each other. 
If we solve the eigenvalue problem of only the first 2x2 
subspace, including the A A and the iVE, we obtain the 
ground state energy, £ = -9.35 MeV (-9.46 MeV), and 
the probability, P NS = 1.57% (2.62%), for the mND s 
(NFs). This clearly means that the couplings between 
the (AA, NE.) and AE components have to be taken into 
account in order to obtain an accurate solution. In the 
case of mNDs, the large coupling potentials, (Vaa-ae) 
and (Vjvh_ae), also enhance the Pne probability. On the 
other hand, for NFs, these coupling potentials hardly en- 
hance Path, though the total binding energy significantly 
increases. The large repulsive energy in the NE compo- 
nent, which mainly comes from the NE—NE diagonal po- 
tential in the 1=0, ^So channel, suppresses the net effect 
of the A A— NE coupling for NFs. On the other hand, the 
NE—NE potential of mNDs is weakly attractive, which 
is consistent with recent experimental data[24j, |25|, |26( . 

In summary, we have performed full-coupled chan- 
nel ab initio calculations for the complete set of doubly 
strange s-shell hypernuclei. Two kinds of YY interac- 
tions, mNDs and NFs, reproduce the APaa^a^ 6 ) 01 
the Nagara event. We obtained bound-state solutions 
for AA H, AA H and A^He by using these YY interactions. 
We thus conclude that a set of phenomenological B S B S 
interactions among the octet baryons in S = 0,-1,-2 
sectors, which is consistent with all of the available ex- 
perimental binding energies of the S = 0, —1, —2 s-shell 
(hyper-) nuclei, can predict a particle stable bound state 
of A aH. For the AA H (and AA He), the probability Path 
by using the mNDs is larger than the Pne by using the 
NF S . We found that the AN - EiV and NE - AE poten- 
tials make a larger Pjve( aa H) for the weak AA — NE and 
attractive NE — NE potentials of mNDs, whereas the net 
effect of the stronger AA — NE coupling potential of NFs 
is suppressed in AA H due to the repulsive NE—NE poten- 
tial in the / — 0, 1 So channel. The one-boson-exchange 
potential models for the PsPs interactions have sound 
bases of the SU (3) symmetry and are widely accepted, 



though uncertainties of the interactions in the S = — 2 
sector are still large because of limitation of experimen- 
tal information. The NSC97 models have a crucial de- 
fect in A A — NE — EE couplings. Therefore, we have 
attempted only two possible cases with the models ND 
and NF, which have different characters in the S = — 2 
sector (strengths of AA — NE coupling and NE — NE 
diagonal potentials). We do hope that a future experi- 
mental facility (e.g., J-PARC) develops our knowledge of 
the S = —2 interactions. 

The authors are thankful to K. Nakazawa for useful 
communications. The calculations were made using the 
RCNP's SX-5 computer and the KEK SR8000 computer. 
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